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ABSTRACT 

We present Chandra and XMM-Newton observations of PLCK G036.7+14.9 from the Chandra-Planck Legacy 
Program. The high resolution X-ray observations reveal two close (~72" = 193 kpc in projection) subclusters, 
G036N and G036S, which were not resolved by previous ROSAT, optical, or recent Planck observations. We 
perform detailed imaging and spectral analyses and use a simplified model to study the kinematics of this system. 
The basic picture is that PLCK G036.7+14.9 is undergoing a major merger (mass ratio close to unity) between 
the two massive subclusters, with the merger largely along the line-of-sight (~80° between the merger axis and 
the plane of the sky from the simplified model) and probably at an early stage (less than ~ 0.4 — 0.7 Gyr since the 
merger began). G036N hosts a small (~27 kpc), moderate cool-core (cooling time f coo i ~ 2.6 — 4.7 Gyr), while 
G036S has at most a very weak cool-core (f coo i ~ 5.7 - 10.3 Gyr) in the central ~40 kpc region. The difference in 
core cooling times is unlikely to be caused by the ongoing merger disrupting a pre-existing cool-core in G036S. 
G036N also hosts an unresolved radio source in the center, which may be heating the gas if the radio source is 
extended. The total mass of the whole cluster determined from XMM-Newton is ~ (5.9 — 8.0) x 10 14 M 0 , and is 
~ (6.7 — 9.9) x 10 14 Mq from Chandra. The Planck derived mass, ~ (5.1 — 6.0) x 10 14 M 0 , is higher than the 
X-ray measured mass of either subcluster, but is lower than the X-ray measured mass of the whole cluster, due to 
the fact that Planck does not resolve PLCK G036.7+14.9 into subclusters and interprets it as a single cluster. This 
mass discrepancy could induce significant bias to the mass function if such previously unresolved systems are 
common in the Planck cluster sample. High resolution X-ray observations are necessary to identify the fraction 
of such systems and correct such a bias for the purpose of precision cosmological studies. 

Subject headings: X-rays: galaxies: clusters —galaxies: clusters: general — galaxies: clusters: individual 
(PLCK G036.7+14.9) — cosmology: observations 


1. INTRODUCTION 

In the concordance Lambda Cold Dark Matter (ACDM) cos¬ 
mology, small scale systems form first and subsequently un¬ 
dergo a series of mergers and accretions to form larger ob¬ 
jects such as groups and clusters of galaxies (e.g., Springel et 
al. 2005; Boylan-Kolchin et al. 2009; Kravtsov & Borgani 
2012). Clusters of galaxies are the largest objects whose in¬ 
ner parts have had time to virialize. Larger systems are still 
forming, with gravity the dominant force governing the behav¬ 
ior of the formation process. On smaller scales, the impact of 
non-gravitational processes such as radiative cooling, feedback 
from supernovae or Active Galactic Nuclei (AGNs) becomes 
significant (e.g., Ponman et al. 1999; Voit et al. 2005). Due to 
their large volumes and high masses, clusters of galaxies can be 
used in a number of ways to constrain cosmological parameters 
(see Allen et al. 2011, and references therein). In particular, the 
number density of clusters as a function of mass and redshift, 
i.e., the mass function, is strongly related to the mean matter 
density Dm and the power spectrum amplitude Cg (e.g., Voit 
2005; Vikhlinin et al. 2009; Mantz et al. 2010). 


The primary link between theory and observation is the mass 
of the cluster. Observationally, the mass can be obtained by 
different techniques with different assumptions and limitations. 
Velocity dispersions of the member galaxies can be used to de¬ 
rive the mass under the assumption of dynamical equilibrium 
(e.g., Binney & Tremaine 1987). However, this method suf¬ 
fers from strong projection effects and is observationally ex¬ 
pensive. In X-rays, by measuring the gas density and temper¬ 
ature profiles, the mass can be derived assuming hydrostatic 
equilibrium (e.g., Sarazin 1988), which can underestimate the 
mass by 10% — 30% due to non-thermal pressure support (e.g., 
Nagai et al. 2007; Piffaretti & Valdarnini 2008; Jeltema et al. 
2008; Battaglia et al. 2012). Strong and weak lensing offer, in 
principle, unbiased mass measurements, but strong lensing al¬ 
lows detailed mass distribution modeling primarily in the cen¬ 
tral regions (e.g., Meneghetti et al. 2013) while weak lensing 
usually has large intrinsic scatter due to its statistical nature 
and the results somewhat depend on the fitting procedure (e.g., 
Meneghetti et al. 2010; Becker & Kravtsov 2011; Rasia et al. 
2012 ). 
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The Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zeldovich 
1972), a distortion of the Cosmic Microwave Background 
(CMB) spectrum produced by inverse Compton scattering 
of the CMB photons as they travel through the intracluster 
medium (ICM), has proven to be very efficient at finding mas¬ 
sive clusters through three recent SZ surveys: the South Pole 
Telescope (SPT; e.g., Reichardt et al. 2013), the Atacama Cos¬ 
mology Telescope (ACT; e.g., Hasselfield et al. 2013), and 
the Planck mission (e.g., Planck Collaboration XXIX 2014). 
Planck is an all-sky survey while SPT and ACT perform deeper 
observations of smaller solid angles with higher angular reso¬ 
lution. In 2011, Planck released an early catalogue of 189 clus¬ 
ter candidates with high reliability based on the first 9 months 
of data (Planck Collaboration VIII 2011). Recently, this early 
sample was increased to 1227 cluster candidates, based on 15.5 
months of data (Planck Collaboration XXIX 2014). The Chan- 
dra-Planck Legacy PrograirQ is collecting Chandra observa¬ 
tions for all the z < 0.35 clusters (165 in total) in the Planck 
early catalogue. Each observation has at least 10,000 source 
counts to ensure sufficient accuracy in characterizing the clus¬ 
ters’ masses, dynamical states, and scaling relations. The ulti¬ 
mate goal is to obtain the local cluster mass function, and com¬ 
pare it to that at higher redshifts (e.g., from SPT and ACT) to 
constrain cosmological parameters. 

In addition, clusters of galaxies are of great interest in an 
astrophysical context. For dynamically relaxed clusters, the 
surface brightness usually exhibits a sharp cusp toward the 
center where the temperature drops below the surrounding re¬ 
gion (e.g., Fabian 1994), a distinct feature of cool-core clus¬ 
ters. For these clusters, whose radiative cooling times in the 
central regions are usually short, some heating mechanisms, 
such as thermal conduction (e.g., Zakamska & Narayan 2003; 
Voigt & Fabian 2004; Sanderson et al. 2009), buoyantly rais¬ 
ing bubbles (e.g., Churazov et al. 2001, 2002; Miraghaei et al. 
2014), shocks (e.g., David et al. 2001; Nulsen et al. 2005), 
sound waves (e.g., Forman et al. 2005; Fabian et al. 2006), 
and cosmic ray leakage (e.g., Guo & Oh 2008; Mathews & 
Brighenti 2008), are believed to operate in order to offset the 
cooling, although the details are poorly known (see McNa¬ 
mara & Nulsen 2007, 2012 for reviews). Cluster outskirts, 
which present opportunities to study a range of physical pro¬ 
cess (see Reiprich et al. 2013 for a review), e.g., deviation 
from hydrostatic/thermal/ionization equilibrium, gas clumping, 
and accretion of the warm-hot intergalactic medium (WHIM), 
are largely unexplored because of their low signal-to-noise ra¬ 
tio. For merging systems, shocks and cold fronts are generated, 
producing contact discontinuities in the surface brightness and 
temperature profiles (e.g., Markevitch et al. 2000; Vikhlinin et 
al. 2001; Markevitch & Vikhlinin 2007). Turbulence is also 
expected, which, together with shocks, might be related to the 
amplification of magnetic fields and acceleration of relativistic 
particles, as inferred from observations of radio halos and radio 
relics (see Feretti et al. 2012, and references therein). Cluster 
merging could also be responsible for the disruption of cool- 
cores (e.g., Sanderson et al. 2006; Burns et al. 2008). Merging 
clusters have also provided direct evidence for dark matter (e.g., 
Clowe et al. 2006) and constraints on the self-interacting dark 
matter cross-section (e.g., Markevitch et al. 2004; Randall et 
al. 2008). 

In this paper, we present Chandra and XMM-Newton ob- 

1 http://hea-www.cfa.harvard.edu/CHANDRA.PLANCK.CLUSTERS/. 

2 http://heasarc.gsfc.nasa.gov/docs/xmin/xmmhp_xmmesas.html. 


servations of PLCK G036.7+14.9 (CIZA J1804.4+1002) from 
the Chandra-Planck Legacy Program. Our X-ray observations 
demonstrate that PLCK G036.7+14.9 is undergoing a major 
merger of two close (in projection) yet separated subclusters. 
PLCK G036.7+14.9 was also observed by ROSAT and in the 
optical band (Ebeling et al. 2002), but was not resolved into 
subclusters by either observation. Due to the size-flux degener¬ 
acy (Planck Collaboration VIII 2011), Planck used the ROSAT 
determined position and size from the Meta-Catalogue of X-ray 
Detected Clusters (MCXC; Piffaretti et al. 2011), so Planck ac¬ 
tually measures the total flux from the whole cluster instead of 
the flux from each individual subcluster. 

This paper is structured as follows: The observations and 
data reduction are presented in Section 2. Section 3 provides the 
imaging analysis, i.e., the morphology and surface brightness 
profile, while detailed spectral analysis, including the treat¬ 
ments of background and other uncertainties, gas mass, total 
mass, and core properties can be found in Section 4. We dis¬ 
cuss our main results in Section 5, with a summary presented in 
Section 6. An Appendix is also supplied for the model 2D/3D 
temperature, gas mass, total mass, and gas mass fraction pro¬ 
files. We adopt the Hubble constant Hq = 10 km s -1 Mpc -1 , 
the matter density parameter = 0.3, and the dark energy 
density parameter = 0.7. At the cluster redshift of 0.1547, 
the luminosity distance is 737.7 Mpc and 1" = 2.682 kpc. Un¬ 
certainties quoted are 1 a throughout this paper. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra Data 

PLCK G036.7+14.9 was observed by Chandra (ObsID 
15098) on 2014 February 5 for 9.6 ks with the Advanced CCD 
Imaging Spectrometer (ACIS-I). The observation was teleme¬ 
tered in VFAINT mode. Standard Chandra data analysis was 
performed with CIAO version 4.6.1 and calibration database 
version 4.6.1. The CIAO tool chandra_repro was applied to 
perform initial processing and to obtain a new event file. Point 
sources were detected by running wavdetect and were excluded 
in further analysis (except for the AGN in the center of the north 
subcluster, see Sections 3 and 4.4). Light curves in the soft- 
and hard-band were examined and no background flares were 
found, so we proceeded with the full exposure. 

The imaging analysis was performed in the 0.7 — 7.0 keV 
band to maximize the signal-to-noise ratio. The period F blank- 
sky background with proper scaling according to the exposure 
time and high energy flux was used. The exposure map was 
weighted by the best fitting thermal model to the spectrum of 
the central part of the cluster. 

In the spectral analysis, the blank-sky background was used 
as the baseline background model. We also varied the blank-sky 
background by ±5% and added a soft-band adjustment to test 
how sensitive the results are to the background (Sections 4.4 
and 4.6). All spectral analysis was performed in the 0.7 — 9.0 
keV band, unless otherwise stated. 

2.2. XMM-Newton Data 

PLCK G036.7+14.9 was observed by XMM-Newton (ObsID 
0692931901) on 2013 March 30 for 10.5 ks (EPIC MOS1), 10.6 
ks (EPIC MOS2), and 9.7 ks (EPIC PN). Full frame mode was 
used for the three cameras with medium filters for the EPIC 


X-ray observations of PLCK G036.7+14.9 


3 


MOS detectors and thin filters for the EPIC PN. Data reduction 
was done with ESAS0 version 5.6 (Snowden & Kuntz 2013), 
which uses SAS version 13.5.0 and calibration files updated to 
2014 February 27. The initial data processing was done by run¬ 
ning emchain and epchain with the default setting in ESAS. 
Out-of-time (OOT) events were also accounted for in the EPIC 
PN analysis. Background flares were identified and rejected by 
applying a 1.5 <7 clipping method to the high energy count rate 
histogram. The resulting “cleaned” exposure times are 6.8 ks 
(EPIC MOS1), 5.4 ks (EPIC MOS2), and 2.0 ks (EPIC PN), re¬ 
spectively. This shows that the XMM-Newton observation was 
affected by high energy particle induced background. Due to 
the limited statistics, we combined data from the three cameras 
in our analysis. Point sources were detected in the 0.5 — 10.0 
keV band and were excluded from further analysis (except for 
the AGN in the center of the north subcluster, see Sections 3 
and 4.4). 

Imaging analysis was performed in the 0.7 —2.0 keV band. 
The count images, quiescent particle background (QPB) im¬ 
ages, residual soft proton (SP) background images (see Section 

4.1.2) , and exposure maps for the three cameras were generated 
and combined with the ESAS task comb. The QPB and SP im¬ 
ages were then subtracted from the combined count image and 
the resulting image was exposure corrected with the combined 
exposure map. 

Spectra from the same regions of the three cameras were 
extracted with ESAS tasks mos-spectra and pn-spectra and 
grouped to contain a minimum of 25 counts per bin. We mod¬ 
eled the various components of the background (see Section 

4.1.2) . Spectra from different regions and from the three cam¬ 
eras were fit simultaneously. The spectral analysis was carried 
out in the 0.7 — 10.0 keV band for EPIC MOS spectra and in 
the 0.7 - 7.0 keV band for EPIC PN spectra. 

3. IMAGING ANALYSIS 
3.1. Morphology 

The XMM-Newton image (Figure 1, left panel) shows that the 
large scale X-ray emission from PLCK G036.7+14.9 is elon¬ 
gated in a northeast-southwest direction. The central region 
reveals two close (^72" = 193 kpc in projection), yet clearly 
separated subclusters (G036N in the north and G036S in the 
south), which suggests that PLCK G036.7+14.9 may be under¬ 
going a merger. The Chandra image (Figure 1, right panel) 
reveals a bow shaped gap (with an angle of ~145°) between 
G036N and G036S, anther indication of interaction between 
the two subclusters. Spectral analysis reveals a hotter region 
between G036N and G036S, confirming that they are interact¬ 
ing (see Section 4.7). Beyond the interaction region, neither 
G036N nor G036S deviates significantly from spherical sym¬ 
metry, suggesting that the merger is probably at an early stage. 

Optical images of PLCK G036.7+14.9 were obtained using 
the MOSaic CAmera (MOSCA) at the 2.56 m Nordic Opti¬ 
cal Telescope (NOT). Using MOSCA in 2 x 2 binned mode 
yielded a pixel scale of 0"217 pixel \ and the field-of-view 
(FOV) is I'.l x 111. Three individual exposures were made 
in each of the Sloan Digital Sky Survey (SDSS) g-, r-, and 
/-band, adding up to total exposure times of 1800 s, 600 s, 
and 600 s, respectively. The images in the r- and /-band were 
smoothed with Gaussian kernels to match the seeing (full width 
at half maximum, FWHM = I "45) of the combined g-band 
image. Object detection and photometry were performed on 

3 For Chandra . the normalized blank-sky background was subtracted; for XMM-i 


the combined image for each filter using the SExtractor soft¬ 
ware (Bertin & Arnouts 1996) in dual image mode with the g- 
band chosen for the reference image. Total magnitudes were 
measured by the MAG_AUTO parameter while colors were 
measured within a 3 1'6 diameter aperture. Photometric zero- 
points were calibrated from SDSS photometry of stars in other 
fields observed using the same setup, and the derived magni¬ 
tudes were corrected for Galactic extinction using the Schlafly 
& Finkbeiner (2011) recalibration of the Schlegel et al. (1998) 
dust maps. The g — r vs. r — i color-color diagram of non- 
stellar objects in the field of PLCK G036.7+14.9 is shown in 
the left panel of Figure 2. A significant overdensity, encom¬ 
passed by the shaded circle, corresponding to the location of 
PLCK G036.7+14.9, is clearly visible. In order to select the 
red sequence of PLCK G036.7+14.9, we fit a linear relation to 
the g — i vs. i color-magnitude diagram (right panel of Figure 2) 
for galaxies within the circle and obtained a best fitting relation 
represented by the dashed line in the right panel of Figure 2. 
The color scatter (1 cr) about this relation, illustrated by the dot¬ 
ted lines in the diagram, is 0.072 mag, which is typical of the 
measured intrinsic scatter around the mean color-magnitude re¬ 
lation for early-type galaxies in rich clusters (e.g., Stanford et 
al. 1998). Galaxies which fall within twice the value of the 
scatter from the best fitting color-magnitude relation, down to 
iAB ~ 21 mag (at fainter magnitudes, the photometric errors be¬ 
come non-negligible compared to the intrinsic scatter around 
the color-magnitude relation), are considered to be the red se¬ 
quence of PLCK G036.7+14.9 and are plotted in red in the 
diagram. The colors of the red sequence are in good agree¬ 
ment with those in Eisenstein et al. (2001) at redshifts close to 
PLCK G036.7+14.9. The Brightest Cluster Galaxies (BCGs) 
of G036N and G036S are located on the red sequence and col¬ 
ored green in particular in both panels of Figure 2. As shown in 
Figure 3 (left panel), the BCGs of the two subclusters, marked 
with blue squares (5" x 5"), are very close to the X-ray centers 
(defined as the X-ray centroids), suggesting that the merger is 
at an early stage and it should be largely along the line-of-sight. 

An X-ray AGN was detected in the 2.0 —7.0 keV band at 
the center of G036N in the Chandra image (Figure 1, right 
panel), while no point source was detected at the center of 
G036S. Interestingly, the NRAO VLA Sky Survey (NVSS) im¬ 
age (right panel of Figure 3) also reveals a radio source, NVSS 
J180431 + 100323, coincident with the center of G036N. This 
radio source should be hosted by the BCG of G036N. The ex¬ 
tent of the radio source is comparable to the angular resolution 
of the NVSS (FWHM = 45"), so it is essentially unresolved. 
The properties of the AGN and the core regions are discussed 
in more detail in Section 4.8. 


3.2. Surface Brightness Profile 

Surface brightness profiles for each subcluster were extracted 
from the background subtracted and exposure corrected im¬ 
ages. We chose a half circular region opposite to the interac¬ 
tion region (see Figure 4) to derive the surface brightness pro¬ 
files for both G036N and G036S. The azimuthally averaged sur¬ 
face brightness profiles are shown in Figure 5. We first tried to 
model these with single j3 models (Cavaliere & Fusco-Femiano 
1976), 

5(r)=5 0 [l + (-) 2 ]- 3/3+1/2 , (1) 

r C 

on, the QPB and SP background images were subtracted (see Section 2). 


4 


Zhang et al. 


where r c is the core radius. We found excess emission in the 
central — 10" regions for both G036N and G036S, so we ex¬ 
cluded the radial bins covering the central — 10" regions and 
refit the surface brightness profiles. The best fitting models 
are plotted in Figure 5, with the best fitting parameters given in 
Table 1. 

As can be seen, ft models describe the surface brightness pro¬ 
files quite well beyond the central — 10" for both G036N and 
G036S, although the reduced X~ for Chandra is larger. The 
best fitting core radii are consistent for Chandra and XMM- 
Newton while the differences between /I values are only 1.9(7 
for G036N and 1.6(7 for G036S. In the central — 10" regions, 
both G036N and G036S show excess emission relative to the ft 
models. Note that the X-ray AGN at the center of G036N is not 
excluded in the surface brightness profiles. Compared to the 
surface brightness profile of G036S, it can be inferred that the 
contribution from the X-ray AGN to the total flux from G036N 
is negligible. We further demonstrate this from the spectral 
analysis in Section 4.4. Given the quality of the data and the 
goodness of the p model beyond the central — 10" regions, we 
did not try more complicated models. 

We also extracted surface brightness profiles in different 
wedge-shaped regions, trying to identify surface brightness dis¬ 
continuities, but the current data do not reveal any unambiguous 
edges. 

4. SPECTRAL ANALYSIS 

4.1. Background Treatment 

4.1.1. Chandra 

The background level is quiescent throughout the Chandra 
observation, so we simply used the period F blank-sky back¬ 
ground as our baseline background model. In addition, to 
check how sensitive the results are to the background variation 
(with time and/or location), we varied the background level by 
±5% (by changing the BACKSCAL keyword by 3=5%), and 
added a soft-band adjustment by fitting the spectrum from a 
(largely) source-free region near the edge of the CCD (Fig¬ 
ure 4, right panel) with an unabsorbed thermal model (best fit¬ 
ting kT = 0.23 keV) allowing the normalization to be negative 
(e.g., Vikhlinin et al. 2005). We treated the uncertainty in the 
background subtraction as a systematic uncertainty (see Section 
4.6). 

4.1.2. XMM-Newton 

As the XMM-Newton data were affected by periods of high 
energy particle induced background (Section 2.2), there may be 
residual SP contamination even after flare screening. Thus we 
decided to model the various components of the background 
as in Snowden & Kuntz (2013). The XMM-Newton back¬ 
ground is very complex (e.g., Nevalainen et al. 2005; Carter 
& Read 2007; Leccardi & Molendi 2008), but can be separated 
into three components: instrumental background, cosmic back¬ 
ground, and solar wind charge exchange (SWCX) background. 

The instrumental background contains QPB (quiescent parti¬ 
cle background), instrumental lines, and residual SP (soft pro¬ 
ton) contamination. The QPB component can be easily sub¬ 
tracted using the filter wheel closed (FWC) data. The instru¬ 
mental lines are dominated by the EPIC MOS Al Ka (— 1.49 
keV), EPIC MOS Si Ka (-1.75 keV), EPIC PN Al Ka (-1.49 
keV), and (six) EPIC PN Cu (—8 keV) lines. We set an upper 
limit of 7.0 keV for EPIC PN spectra, leaving only three instru¬ 
mental lines, which were modeled as three Gaussians with zero 


width. The residual SP contamination was modeled as a power 
law not folded through the instrumental effective areas, with the 
indices linked together for the EPIC MOS cameras. 

The cosmic background was modeled as three components: 
an unabsorbed thermal component with kT = 0.1 keV, repre¬ 
senting the emission from the Local Hot Bubble (LHB) or he¬ 
liosphere; an absorbed thermal component, allowing the tem¬ 
perature to vary, representing the emission from the hotter 
Galactic halo and/or the intergalactic medium; and an absorbed 
power-law with the index fixed at 1.46 representing the back¬ 
ground from unresolved cosmological sources. Both the ab¬ 
sorbed and unabsorbed thermal components have zero redshift 
and solar abundance. 

The SWCX background only produces line emission, with 
the strongest lines being O VII (—0.57 keV) and O VIII (—0.65 
keV) lines. We set a lower energy cutoff of 0.7 keV, so that the 
SWCX background can be neglected. 

The normalizations of the various components of the back¬ 
ground were determined by simultaneously fitting the spec¬ 
tra from the global cluster regions (see Figure 4). We define 
seven regions for the global spectral fitting: core (small circle, 
35" in radius), cluster (ellipse), cluster minus core, for both 
G036N and G036S, and a large annulus (480"-720" in radii) 
free of cluster emission. The normalizations of the cosmic 
background components, together with the normalization of the 
source spectrum (a single absorbed thermal model), were com¬ 
puted in units of arcmin 2 . All the normalizations of the back¬ 
ground components were left as free parameters in the global 
spectral fitting, and were fixed in the subsequent spectral anal¬ 
ysis (e.g., temperature profile), accounting for the solid angle. 
The temperature of the absorbed cosmic background compo¬ 
nent and the indices of the residual SP contamination were also 
fixed at the best fitting values in the subsequent spectral analy¬ 
sis. 

4.2. Hydrogen Column Density 

Since PLCK G036.7+14.9 is located at low galactic latitude 
(b = 14.9°), we expect the total Hydrogen column density, Nh , 
to be higher than the Leiden/Argentine/Bonn (LAB) H I sur¬ 
vey (Kalberla et al. 2005) value of Vhi = 0.087 x 10 22 cm~ 2 , 
due to additional absorption by molecular Hydrogen. Indeed, 
the best fitting Nh obtained by the global spectral fitting is 
0.092+QQQg x 10 22 cm' 2 for XMM-Newton, 0.164+°;°[ 2 x 10 22 
cm- 2 for Chandra (baseline background model and baseline 
background model varied by ±5%), and 0.195town x 10 22 
cm- 2 if a soft-band adjustment is added to the Chandra base¬ 
line background model. 

XMM-Newton and Chandra give significantly different Nh, 
by almost a factor of two. To investigate the origin of this dis¬ 
crepancy, we refit the spectra with energy cutoffs of 0.3 keV, 
0.5 keV, 0.7 keV and 1.5 keV. The best fitting Nh values are all 
consistent for each observatory, but still significantly different 
between the two observatories. We also tried different abun¬ 
dance tables, and found that while Nh changes with different 
abundance tables, the ratio between the two observatories re¬ 
mains almost unchanged. Therefore, we conclude that the dis¬ 
crepancy is most likely caused by cross-calibration issues (see 
Section 4.6 on their effects on the total mass determinations). 

To examine whether Nh is different between G036N and 
G036S, and whether Nh changes with radius, we fit all spec¬ 
tra (Section 4.4) of G036N and G036S with Nh free in each 
annulus. We found that all the Nh values are consistent within 
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the lcr errors for Chandra, and consistent within the 2cr errors 
for XMM-Newton. Therefore, in our spectral analysis we linked 
Nh in different regions together. 

In addition, we estimated the total Nh with the relatiorQ 

A^H.tot = Mil + 2A^H 2 = Nm + 2A^h 2 ,max (1 — e Wc ) , (2) 

where .max = 7.3 x 10 20 cm“ 2 , E(B — V) = 0.162 is the dust 
extinction, N c = 3.0x 10 2() cm 2 , and a = 1.1 (Willingale et al. 
2013). We obtained Ah, tot = 0.136 x 10 22 cm -2 . All the results 
obtained with Nh fixed at this total value are also included for 
comparison. 

To summarize, we used Nh = 0.092 x 10 22 cm~ 2 {or XMM- 
Newton, Nh = 0.164 x 10 22 cm -2 for the Chandra baseline 
background model and the baseline background model varied 
by ±5%, and Nh = 0.195 x 10 22 cm~ 2 for the Chandra base¬ 
line background model with a soft-band adjustment. We also 
included Nh = 0.136 x 10 22 cm~ 2 for both observatories for 
comparison. See Table 2 for the Nh values used for the two ob¬ 
servatories and for the different Chandra background models. 
We treated the uncertainty in Nh as a systematic uncertainty 
(see Section 4.6). 


4.3. Redshift 

The only redshift available for PLCK G036.7+14.9 is from 
the Cluster in the Zone of Avoidance (CIZA) project, which 
measured a spectroscopic redshift z = 0.1525 (Ebeling et al. 
2002). However, PLCK G036.7+14.9 was treated as a single 
cluster in Ebeling et al. (2002), so we decided to use our best 
fitting X-ray redshifts for G036N and G036S. 

Table 2 lists the redshifts obtained by XMM-Newton and 
Chandra with different Chandra background models and dif¬ 
ferent Nh values. As can be seen, a different Nh for the same 
observatory and the same background model (Chandra) only 
has a minor effect on the derived z . The obtained best fitting 
redshifts differ by <2% and are consistent. However, XMM- 
Newton and Chandra give values that differ by ^13% (^3cr) 
and ~7 % (~2cr) for G036N and G036S, respectively, which 
are probably caused by cross-calibration issues between the two 
observatories. 

We simply averaged the best fitting redshifts over XMM- 
Newton and Chandra (baseline background model), and those 
obtained with different Nh values, to obtain z n = 0.1501 ± 
0.0022 and z s = 0.1592 ±0.0020, where (and hereafter) we use 
the subscripts n and s to distinguish between the north and south 
subclusters. We assume that the difference, which has a signif¬ 
icance of 3cr, between the redshifts of the two subclusters is 
caused by their relative motion, and the redshift of the cluster 
is the mean of the two, 0.1547 ±0.0015, which is very close to 
the optical redshift of 0.1525. All the distance-relevant quanti¬ 
ties are relative to a redshift of 0.1547. In our spectral analysis, 
the redshift is fixed at the value shown in Table 2 according 
to the context. The uncertainty in the redshift was treated as a 
systematic uncertainty (see Section 4.6). 

4.4. Temperature Profile and Deprojection 


Azimuthally averaged temperature profiles were derived in 
half circular regions opposite to the interaction region (see Fig¬ 
ure 4) for G036N and G036S. Spectra were extracted in the 
0-35" (0-35"), 35"-70" (35"-70"), 70"-130" (70"-140"), 130"- 
260" (140"-280"), 260"-600" (280"-600") regions for G036N 
(G036S) for XMM-Newton, and in the 0-35'Q (0-45"), 35"- 
80" (45"-100"), 80"-200" (100"-300"), 200"-450" regions for 
G036N (G036S) for Chandra. All spectra for each observatory 
were fit simultaneously. Due to the limited statistics, the abun¬ 
dances in each annulus cannot be constrained, so we fixed them 
at the global best fitting values. 

The temperature profiles for each subcluster and for each set 
of {Nh, z} values in Table 2, are displayed in Figure 6. In each 
panel of Figure 6, the temperature profiles for XMM-Newton 
and for Chandra with different background models (see Section 
4.1.1) are indicated by different symbols. When using the same 
Nh value, XMM-Newton and Chandra generally give consistent 
results, although at large radii (> 200"), a direct comparison is 
not straightforward as they do not probe exactly the same re¬ 
gion. For G036N, there is an indication of a temperature drop in 
the central region, although the uncertainties are large. At large 
radii, the temperature profile declines, dropping to ~\/2 (~l/4) 
of the “peak” value for Chandra ( XMM-Newton ). This behavior 
is qualitatively consistent with large samples of cluster temper¬ 
ature profiles (e.g., Vikhlinin et al. 2005; Pratt et al. 2007; Lec- 
cardi & Molendi 2008). For G036S, XMM-Newton exhibits a 
temperature drop in the innermost bin, and a declining temper¬ 
ature profile at large radii, while Chandra, in the much smaller 
radial range it covers because the cluster is close to the chip 
boundary, gives an essentially isothermal temperature profile, 
but with large uncertainties. Different Nh values and different 
Chandra background models produce the same trend as previ¬ 
ously stated for the temperature profiles for both G036N and 
G036S. 


The observationally determined temperature profile is the 
projected 2D profile, which contains emission from different 
parts of the cluster along the line-of-sight. We used the func¬ 
tional form, 


73dM = To 


( _ l _j^cool . 

V ''cool > 


J min 

To 


(~y a 

y r, > 




(3) 


proposed by Vikhlinin et al. (2006), to model the 3D temper¬ 
ature profile. Assuming spherical symmetry, this 3D tempera¬ 
ture profile is then projected along the line-of-sight, weighted 
by n 2 /T 3 / 4 5 6 (Mazzotta et al. 2004), where n e is the electron den¬ 
sity and T is the 3D temperature, and is fit to the observed 2D 
temperature profile. As this functional form has many free pa¬ 
rameters, it can describe the temperature profile very well. The 
results are shown in the Appendix. 


4.5. Gas Mass, Total Mass, and Gas Mass Fraction 

For a surface brightness profile described by a model, the 
electron density is given by 

n e {r) =«<.,o[l ±(—) 2 r 3/3/2 . (4) 

r c 

The central dcnsityQ, n e{] , was determined by fitting the spec¬ 
trum extracted in a 200" circular region, utilizing the relation 


4 http://www. swift.ac.uk/analysis/nhtot/index.php. 

5 We did not exclude the AGN at the center of G036N, since the bolometric X-ray luminosity (0.01 — 100 keV in practice) of the AGN is only ~1% compared to that 
of the first bin, so the presence of the AGN does not affect the temperature of the first bin for G036N. 

6 Note that there is excess emission relative to the /3 model (see Figure 5) in the central ~10 // region, so the “central density”, n e o, is not equal to the physical central 
density. We use “central density” just in terms of the functional form. However, it is still accurate to use n e $ to determine the gas mass at large radius, as the central 
~10 // = 27 kpc region only occupies a tiny volume. 
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between the normalization K of the thermal model and electron 
number density n e , 

1(T 14 f 

K= ^m Tijp () 

where Da is the angular diameter distance, nn is the Hydrogen 
number density, and the integration is over a cylindrical vol¬ 
ume. The gas mass within radius r is then 

M ? (r) = [ 4KHen e flm p [l + (—) 2 ]~ 3/3/2 r 2 c/r, (6) 

Jo r c 

where fi e = 1.172 is the mean molecular weight of electrons, 
and m p is the mass of a proton. 

The total mass within radius r can be derived, under the as¬ 
sumption of hydrostatic equilibrium, from (e.g., Sarazin 1988) 


M(r) 


kTr 
G finip 


, d In n„ 

( - 1 

y d In r 


d\nT 

dlnr 


(7) 


where k is the Boltzmann constant, G is the gravitational con¬ 
stant, p = 0.614 is the mean molecular weight of the total par¬ 
ticles, and T is the 3D temperature at radius r. Since PLCK 
G036.7+14.9 is undergoing a merger (Sections 3.1 and 4.7), we 
avoided the interaction region and measured the surface bright¬ 
ness and temperature profiles in the outer halves of both sub¬ 
clusters (see Figure 4). The surface brightness profiles from 
these regions are well described by j3 models (Figure 5 and Ta¬ 
ble 1), and the morphologies do not deviate significantly from 
spherical symmetry (Figure 1), so hydrostatic equilibrium in 
these regions should be a reasonable approximation. 

Uncertainties were estimated from Monte Carlo simulations. 
The errors of the measured 2D temperature profile, the param¬ 
eters of the surface brightness profile, and the normalization 
K of the spectral fit to the 200" region, are assumed to obey 
Gaussian distributions. For each Monte Carlo simulation, a set 
of randomly drawn 2D temperature profile, jj model parame- 
terfl and normalization K from their respective Gaussian dis¬ 
tributions, was used to obtain the gas mass, total mass, and gas 
mass fraction profiles, using the deprojection method in Sec¬ 
tion 4.4. This process was repeated 1000 times. In practice, too 
steep temperature profiles often lead to unphysical mass deter¬ 
minations. We only accept physical solutions with p tot > p gas , 
where p tot and p gas are the total and gas densities respectively, 
and which are convectively stable, d In T/d In p gas <2/3. Fi¬ 
nally, we quote the mean of all accepted solutions as the “best 
fitting” gas mass, total mass, and gas mass fraction, with the 
uncertainties as the 1 a standard deviation of all the accepted 
solutions. 

In particular, we are interested in the gas mass, total mass, 
and gas mass fraction inside three characteristic radii, r\ (A = 
2500,500,200), where r& is the radius within which the mean 
matter density of the Universe is A times the critical density 
at the cluster redshift. The results are listed in Table 3. Since 
G036S is located near the chip boundary in the Chandra FOV 
(Figure 1), and the Chandra temperature profile does not cover 
the declining region at large radii (Figure 6), we did not mea¬ 
sure the mass-related quantities for G036S with Chandra. 


4.6. Systematic Uncertainties 


within their lcr errors for some parameters which affect the 
mass determinations. We treat the uncertainties in Nh, z, 
other cross-calibration uncertainties between XMM-Newton 
and Chandra, and the uncertainties in Chandra background 
subtraction, as systematic uncertainties. For each set of {Nh, 
z} in Table 2, new temperature profiles for G036N and G036S 
were derived. The characteristic radii, gas mass, total mass, and 
gas mass fraction were then calculated following the method 
given in Section 4.5. We also tried different background sub¬ 
traction methods for the Chandra data (see Section 4.1.1). The 
results are presented in Table 3, from which we can draw the 
following conclusions: 

1. Comparing Nh fixed at the total value to the best fitting 
values, all the quantities differ by less than 1.6(7. The differ¬ 
ences between the “best fitting” values of r&, M g ^, M\, and /a 
for XMM-Newton are < 8%, < 7%, < 25%, and < 16%, while 
those for Chandra are < 10%, < 8%, < 31%, and < 27%, re¬ 
spectively. The larger differences for Chandra are caused by 
the Chandra baseline background model with a soft-band ad¬ 
justment, without which the differences for all four quantities 
for Chandra are reduced by a factor of at least 2. 

2. Comparing XMM - Ne w ton to Chandra, all quantities differ 
by no more than 2.5a. Varying the Chandra baseline back¬ 
ground by +5% (back*0.95) gives the closest results to XMM- 
Newton, with < 9%, < 6%, < 28%, and < 20% differences 
between the “best fitting” values of r&, M gj a, Ma, and /a, 
while adding a soft-band adjustment to the Chandra baseline 
background gives the largest differences to XMM-Newton, with 
< 34%, < 27%, < 154%, and < 46% differences for the above 
four quantities, respectively. 

3. Comparing Chandra different background models, all the 
quantities are consistent within their la errors. The differences 
(with respect to the baseline background model) of r\, M g & , 
Ma, and /a when varying the baseline background model by 
±5% are < 6%, < 5%, < 15%, and < 10%, while adding a 
soft-band adjustment to the Chandra baseline background gives 
much larger differences, < 19%, < 19%, < 75%, and < 28%, 
respectively. 

By comparing the results in the above mentioned cases, we 
can see that adding a soft-band adjustment to the Chandra base¬ 
line background always produces the largest differences of the 
“best fitting” values of r a, M gj a, Ma, and /a. Actually, adding 
a soft-band adjustment also increases the best fitting Nh for 
the global cluster (Table 2), which further affects the temper¬ 
ature measurements. This “soft-band-adjustment-V/r-k7”’ de¬ 
pendence makes the mass determination much more uncertain, 
which gives rise to the largest differences of the “best fitting” 
values for the mass-related quantities in Table 3. Without accu¬ 
rate Nh determination, we consider this model less constraining 
and do not include it in our further analysis and final results. 

To summarize, uncertainties in Nh, Z, other XMM-Newton 
and Chandra cross-calibration uncertainties, and Chandra 
background subtraction, do not lead to any significantly dif¬ 
ferent (< 2.5a) mass determinations. XMM-Newton gives a 
mass M 200 ~ (3.1 —4.2) x 10 14 M 0 for G036N and M 200 ~ 
(3.0 —4.7) x 1O 14 M 0 for G036S, while Chandra gives a some¬ 
what higher mass, M 200 ~ (3.4 —6.7) x 10 14 M Q , for G036N. 


From previous sections, we see that, due to cross-calibration 
issues, XMM-Newton and Chandra give results inconsistent 

7 It is worth noting that, strictly speaking, r c and j3 are not independent (e.g.. 
overestimated. 


4.7. Is There Shocked Gas between the Subclusters? 
Henriksen & Tittley 2002), so the errors of the gas mass obtained in this way are 
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The X-ray morphology of PLCK G036.7+14.9 suggests that 
it is undergoing a merger between G036N and G036S (Section 
3.1). To confirm this spectroscopically, we extracted spectra 
from a 130 ,, x35 ,, rectangular region (Region 1, see Figure 7) 
between G036N and G036S, and from the symmetric regions 
with respect to their centers (Region 2 for G036N and Region 
3 for G036S), and fit them to single thermal models. We found 
that the temperature in Region 1 is indeed higher (by a factor of 
~2 based on the Chandra data) than those in Regions 2 and 3, 
which could be interpreted as shock heating during the merger. 

Using the Rankine-Hugoniot jump conditions, the Mach 
number, is related to the temperature jump as (e.g., Sarazin 
2002) 

1 _ 2 7-1 
C ^ 2 (y+1) + y+ 1 ’ 

and 

i = [I(Z±l)2(Z2_ 1 )2 + Zk]l/2_IZ±l(Zk_ 1 ) (9) 

c 4 y — D 7j ’ 7i 2 y— 1 7) 
where the subscripts 1 and 2 denote the preshock and postshock 
gas, C = P 2 /P 1 is the shock compression, and y = 5/3 is the 
adiabatic index. If we simply assume that the temperature of 
the preshock gas is given by the gas temperature in Region 2 
(Region 3) and that of the postshock gas is given by the gas tem¬ 
perature in Region 1, using the above relations, we can estimate 
the Mach number of the shock propagating in G036N (G036S) 
both for XMM-Newton and Chandra (see Table 4). Note that 
due to projection effects, the inferred Mach number is only a 
lower limit. 

XMM-Newton suggests a weak shock ~ 1.0 — 1.6), while 
Chandra indicates a stronger shock (.X/ ~ 1.5 — 2.7). The dif¬ 
ference is caused by, when comparing the same Nh, the lower 
temperature in the interaction region measured from XMM- 
Newton (7.65 — 9.19) keV compared to Chandra (10.52— 16.86 
keV). Actually, from the Highest X-ray FLUx Galaxy Cluster 
Sample (HIFLUGCS; Reiprich & Bohringer 2002), Schellen- 
berger et al. (2014) found systematic lower temperatures given 
by XMM-Newton compared to Chandra for high temperature 
systems, especially for >8 keV clusters. When comparing the 
Mach number obtained with the best fitting Nh values, the large 
difference in Mach number between XMM-Newton and Chan¬ 
dra is then caused by higher temperatures in the preshock re¬ 
gions (Regions 2 and 3) given by XMM-Newton. This is be¬ 
cause the best fitting XMM-Newton Nh is almost a factor of 2 
less than that of Chandra (Section 4.2), and a smaller Nh in¬ 
creases the temperature, leading to smaller temperature differ¬ 
ence in the postshock region, and larger temperature difference 
in the preshock regions. Thus, cross-calibration issues again 
manifest themselves in terms of the Mach number discrepancy. 


4.8. Subcluster Cores 


From Section 3.1, we know that G036N has an X-ray AGN in 
the center. We extracted a spectrum in a circular region with a 
radius of 3", and used the surrounding region (4" — 8" annulus) 
as the background. By fixing the index of the power-law model 
at 1.7, we acquired a bolometric X-ray luminosity (0.01 — 100 
keV in practice) of (1.45 — 1.50) x 10 42 erg s _1 . 

The surface brightness profiles (Figure 5) and the tempera¬ 
ture profiles (Figure 6) suggest that the gas may be cooling in 
the central regions of both G036N and G036S. The isobaric 
cooling time was calculated as 


H 

Loot — 7 

^bol 


5 

2 M'» P 
f-'bol 


( 10 ) 


where H is the enthalpy, L ho i is the bolometric X-ray luminos¬ 
ity, and M gas is the gas mass. We then compared this to the look 
back time to z = 1, which is tr = 7.7 Gyr, a representative time 
for a cluster to relax and develop a cool-core (e.g., Birzan et al. 
2004; Hudson et al. 2010). If f coo i < 0., the gas may be cooling. 
Due to the much larger PSF of XMM-Newton , we only used the 
Chandra data to study the cores. From Figure 4, varying the 
blank-sky background by ±5% has little effect on the tempera¬ 
ture in the central bright part of the subclusters, so we only used 
the baseline background model to study the cores. 

For G036N, the cooling time (single thermal model) in the 
central 3" — 10" (full annulus) region is f C ooi = 2.64 — 4.71 Gyr, 
while that in the outer 10" — 35" (half annulus to avoid the in¬ 
teraction region) region is 6.82 — 12.15 Gyr, so only the gas in 
the central ~10" (27 kpc) region may be cooling. The abun¬ 
dance obtained with a single thermal model is 0.85 — 2.44 Zq, 
significantly larger than the global abundance (0.09 — 0.47 Z Q ) 
of G036N, suggesting enrichment from the BCG. A two ther¬ 
mal components model, with the higher temperature compo¬ 
nent fixed at the cluster temperature and abundance, does not 
improve the fit and results in essentially zero normalization 
for the higher temperature component. Adding a cooling flow 
model to a single thermal model does not improve the fit either, 
and the mass deposition rate given by the normalization of the 
cooling flow model is highly uncertain with M spe c = 2.3 —27.6 
M 0 yr^ 1 . 

For G036S, the cooling time (single thermal model) in the 
central 0" — 15" (full annulus) region is t coo i = 5.65 — 10.27 
Gyr, already comparable to or greater than f/ , so cooling should 
be very weak in this region. Adding a cooling flow model, we 
obtained a lcr upper limit on the mass deposition rate of 14.3 
M© yr^ 1 . 


5. DISCUSSION 
5.1. The Merging Activity 

Previous ROSAT and optical studies (Ebeling et al. 2002) did 
not resolve the morphological details of PLCK G036.7+14.9. 
With our high resolution Chandra and XMM-Newton obser¬ 
vations, two close (in projection) subclusters, G036N and 
G036S, were clearly resolved, and spatially resolved spec¬ 
troscopy could be performed. Spectral analysis produces very 
similar redshifts for G036N and G036S (Table 2), and reveals a 
higher temperature in the region between them compared to the 
symmetric regions away from their centers which could be due 
to shock heating during the merger (Section 4.7). 

The morphologies of G036N and G036S excluding the inter¬ 
action region do not deviate significantly from spherical sym¬ 
metry (see Figure 1) and their surface brightness profiles are 
well described by single p models beyond the central ~10" 
regions (see Figure 5), suggesting that the merger is probably 
at an early stage so that perturbations generated by the merger 
have not yet disturbed the outer regions. Another piece of ev¬ 
idence for an early stage of the merger comes from the good 
spatial correspondence between the X-ray emission centroids 
of the two subclusters and their BCGs (see left panel of Figure 
3). This also favors a merger close to the line-of-sight. 

Projection effects hamper our knowledge of the geometry of 
the merger. The only observables are the projected distance, 
d p = 193 kpc, and the radial velocity difference, v r = 2356 km 
s~ 1 , between the two subclusters. To gain some insight into the 
merger kinematics, we applied the simplified model in Ricker & 
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Sarazin (2001) (see also Sarazin 2002; Sauvageot et al. 2005) 
to G036N and G036S. This model treats the two subclusters 
as point masses, and assumes that they initially expand away 
from each other in the Hubble flow. Due to their gravitational 
attraction, they collapse after reaching the greatest separation 
do, and finally merge. Let Coal be the age of the universe when 
core coalesce occurs. The parameters do and Coal are related by 
Kepler’s Third Law as 

4~[2G(M„+M i .)] 1 /3(^)2/3 ) (11) 

7T 

where M is the total mass (assumed to be M 200 , see Table 3). Let 
d be the 3D distance between the two subclusters. The merger 
velocity v is then 

v 2 « 2 G(M n +M,)(1 - ^)[1 - (12) 

where b is the impact parameter as defined in Sarazin (2002), 
and is assumed to be zero in our case since it is usually very 
small compared to do (Sarazin 2002). The parameters d and v 
are related to the observables by 


d= dp 
u n •) 

cos 6 

(13) 

V r 

' “ sin 9 ’ 

(14) 


where 9 is the angle between the merger axis and the plane of 
the sky. On the other hand, since the two subclusters have not 
attained core coalesce, Coal can be written as 


Coal ' 


' Cge T ' 


(15) 


where Cge = 11-52 Gyr is the age of the Universe at the clus¬ 
ter redshift in our adopted cosmology. To obtain the time since 
the merger began, f merge , he., the time from when the two sub¬ 
clusters started physical contact to the present day, we have to 
treat the two subclusters as extended sources and f mer ge can be 
approximated by 


t merge 



(16) 


where r is the virial radius (assumed to be noo, see Table 3). 
Solving Equations (11) — (16), we obtained 9 ~ 80°, t meTge ~ 
0.8 Gyr (see Table 4). Although overly simplified, this model 
produces consistent results with our expectations, i.e., an early 
stage of the merger, mostly along the line-of-sight, which can 
explain the main observational features outlined in the above 
two paragraphs. Targeted simulations are required to study in 
detail the dynamics of this system, which is beyond the scope 
of this paper. 

To summarize, PLCK G036.7+14.9 is undergoing a major 
merger by two nearly equal mass (~ 1 : 1) subclusters. The 
merger should be at an early stage, and should be happening 
mostly along the line-of-sight. A simplified model suggests 
that the merger probably began ~0.8 Gyr ago, with an angle 
between the merger axis and the plane of the sky of ~80°. 


5.2. X-ray Derived Mass vs. SZ Derived Mass 

In Sections 4.5 and 4.6, we derived the total masses of 
G036N and G036S under the assumption of hydrostatic equi¬ 
librium of the ICM using the regions not affected by the on¬ 
going merger. Chandra does not cover a large enough region 


to measure the total mass of G036S. If we assume that Chan¬ 
dra would give the same mass ratio of G036N to G036S as 
XMM-Newton, and the relative errors are the same for G036N 
and G036S for Chandra, we can obtain the total mass of 
G036S that would be “measured” by Chandra. Adding the to¬ 
tal X-ray masses of G036N and G036S together, we obtained 
the total X-ray derived mass (measured inside the X-ray de¬ 
rived rsoo, Of 500 , see Table 3) of PLCK G036.7+14.9, which 
is Mx _500 = (5.91 —8.00) x 10 14 M 0 for XMM-Newton and 
Mx ,500 = (6.66 — 9.85) x 10 14 M 0 for Chandr$\. 

The integrated SZ signal, Ysz, is proportional to the total elec¬ 
tron pressure of the cluster, 

D\Y SZ = f Pe dV, (17) 

m e c z J 

where Or is the Thomson cross-section, m e is the mass of an 
electron, c is the speed of light, and p e is the electron pres¬ 
sure. Due to the size-flux degeneracy (Planck Collaboration 
VIII 2011), additional information is needed to refine the blind 
Planck SZ flux measurement. Using the X-ray measured cluster 
position and size can break the degeneracy, and if the redshift 
of the cluster is known, it may be used to further break the de¬ 
generacy (Planck Collaboration XXIX 2014). After obtaining 
the refined SZ flux, the total SZ mass, M S z, can be estimated 
using the Malmquist-bias corrected scaling relation between 
Ygz and M sz (Planck Collaboration XX 2014). With the above 
method, Planck Collaboration XXIX (2014) published the SZ 
flux, Tsz, 500 , and SZ mass, Msz, 500 , measured inside the SZ 
rsoo, rsz.soo, for 813 confirmed clusters with measured redshifts 
out of the 1227 all-sky SZ sources using the first 15.5 months 
of data. For PLCK G036.7+14.9, Y sz ,500 = 0.00234+°;““? 
arcmin 2 (adopting aredshift of 0.1525) and Msz ,500 = 5.54+Q+ 2 
(5.11-5.96) xlO 14 M 0 . 

It is clear that Msz ,500 is greater than the X-ray mass of ei¬ 
ther G036N or G036S (see Table 3), which is fully expected 
since Planck did not resolve the two subclusters so the derived 
SZ mass is the mass of the whole cluster. However, compared 
to the X-ray mass of the whole cluster, Mx. 500 * especially the 
Chandra “measurement” (note that Chandra actually does not 
directly measure the total mass of the whole cluster), Msz ,500 is 
smaller, although there is a small overlapping range between 
Msz ,500 and the XMM-Newton derived mass because of the 
large uncertainties in the X-ray mass measurements (however, 
see Footnote 8). 

In the derivation of the SZ mass, Planck Collaboration XXIX 
(2014) used the best fitting Tjz.soo — Msz ,500 relation of Planck 
Collaboration XX (2014). Even if we accounted for the uncer¬ 
tainties in that relation, and chose the parameters that maximize 
or minimize Msz, 500 , we obtained the most conservative SZ 
mass estimate, Msz,sot) — (4.99 — 6.06) x 10 14 M 0 . This only 
slightly increases the overlapping range between the SZ and 
XMM-Newton mass measurements, while the lower SZ mass 
compared to X-ray mass determinations remains unexplained. 

The difference between the X-ray and SZ masses should 
not be entirely caused by the different extraction regions for 
the SZ and X-ray mass measurements. In Figure 8, we show 
rsz .500 used by Planck (Planck Collaboration XXIX 2014; Pif- 
faretti et al. 2011; Ebeling et al. 2002) and rx.soo used by 
XMM-Newton and Chandra. The SZ extraction region is larger 
(by ~44%) than the XMM-Newton extraction regions, but the 


8 Note that the values quoted are simply the lowest and highest masses for different Nh and different background models ( Chandra ) in Table 3, so the uncertainties 
are the most conservative estimates. This is important when comparing to the SZ derived mass. Also note that the large uncertainties in the X-ray mass are dominated 
by systematic uncertainties (see Table 3). 
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XMM-Newton measured mass is higher. The Chandra extrac¬ 
tion regions are close to (~3% smaller than) the SZ extraction 
region, but the Chandra measured mass is even higher. There¬ 
fore, while different extraction region is a factor affecting the 
mass determinations from X-ray and SZ measurements, some 
other factors may be more important to explain the difference. 

As we show in Sections 5.1 and 5.3, the merger is most likely 
at an early stage. The SZ signal Y$z relates to the mass M 
via a power-law, which can be written as Y$z = AM^, where 
A and rj can be considered as constants for our purpose. For 
an early stage of the merger, the total mass of the whole cluster 
is M„ +M S . However, if the two subclusters are not resolved, 
as in the case of the Planck observation, the mass of the whole 
cluster determined in this way will be (M,7 +M i I, ) 1//r| . From 
Table 3, M„ ss M v ; rj = 1.79 as used by Planck Collaboration 
XXIX (2014). This gives a total mass a factor of ^0.74 smaller 
than M„ +M S , which is about the difference between the SZ 
derived mass and the X-ray derived mass. This demonstrates 
that the mass discrepancy between the SZ and X-ray measure¬ 
ments is caused by the fact that Planck does not resolve the two 
subclusters and interprets the whole system as a single cluster. 

Based on simulations, Kravtsov et al. (2006) propose that 
the X-ray analogue of Ysz, Yx = M gas x Tx, where Tx is the 
emission weighted temperature within a certain region, is a low 
scatter mass proxy and less sensitive to the cluster dynamical 
state. We computed the total Tv.soo by adding those of G036N 
and G036S, where Tx is measured in the [0.15 — 0.75]rsoo re¬ 
gion and M gas is measured within r $oo (see Table 5 for details). 
We list Tx,500 in Table 5. Tx.soo and Mx.soo follow the scaling 
relations presented in Arnaud et al. (2010). 

Neglecting gas clumping, Yx relates to Ysz by (e.g., Arnaud et 
al. 2010; Rozo et al. 2012; Planck Collaboration XXIX 2014) 

D\Ysz _ QT 1 <n e T> _ c <n e T> 

Y x m e c 2 p e m p <n e >T x <n e >T x 

Therefore, the ratio D 2 A Ysz/CYx should be close to 1, although 
the exact value depends on the structure of the cluster. Here¬ 
after, we will refer Ysz/Yx to D\Y$z/CYx for simplicity. Previ¬ 
ous studies show that Tsz,50o/Tx,500 is between 0.8 and 1 (e.g., 
Arnaud et al. 2010; Rozo et al. 2012; Planck Collaboration 
XXIX 2014). Our results are consistent with these studies, al¬ 
though with large uncertainties (see Table 5). 

In summary, the SZ mass given by Planck is higher than the 
X-ray derived mass by XMM-Newton and Chandra for either 
subcluster, but is lower than the X-ray mass of the whole clus¬ 
ter, due to the fact that Planck does not resolve the two sub¬ 
clusters and interprets the whole system as a single cluster. The 
Ysz/Yx ratio is consistent with previous studies, although with 
large errors. 

5.3. Cool-Core, Merger, and Feedback 

The surface brightness profiles show excess emission in 
the central ~ 10 ” regions relative to single p models for both 
G036N and G036S (see Figure5), and the temperature profiles 
show indications of temperature drops toward their centers (see 
Figure 6 ). The cooling time of the central 3” — 10” region for 
G036N is 2.64 — 4.71 Gyr, and beyond this region, the cooling 
time is longer than 7.7 Gyr. The cooling time of the 0 — 15” 
region for G036S is 5.65 — 10.27 Gyr, comparable to or longer 
than 7.7 Gyr (see Section 4.8 for details). Therefore, G036N 
has a small (~10” = 27 kpc), moderate cool-core, while G036S 


has at most a very weak cool-core within the central ~15” = 40 
kpc region. The spectroscopically derived mass deposition rate 
(see Section 4.8) confirms this scenario. 

It is interesting to ask what causes the difference between the 
cores of the two subclusters. One possibility is that the ongoing 
merger has disrupted a pre-existing cool-core in G036S (e.g., 
Sanderson et al. 2006; Burns et al. 2008; Rossetti & Molendi 
2010; but see, e.g., Poole et al. 2008). From Table 4, the merger 
began ~ 0.78 — 0.80 Gyr ago based on the simplified merger 
model (Section 5.1), while the time for the shock to propagate 
to the center of G036N (G036S), t r _„ (t r _ s ), t r = r 2 oo/-^c SOU nd> 
is 0.65 — 1.12 (0.67 — 1.05) Gyr obtained from XMM-Newton, 
and is 0.43 — 0.74 (0.41 — 0.67) Gyr from Chandra, respec¬ 
tively. Taken at face value, Chandra favors a scenario in which 
the shock has passed the cores of the two subclusters, while the 
situation for XMM-Newton is less clear. 

To further examine this question, we divided the central 
3" — 10” annular region of G036N into two equal area half an¬ 
nular regions, with the symmetry axis perpendicular to the line 
connecting the centers of the two subclusters. We compared the 
derived temperatures in the two half annular and full annular re¬ 
gions, finding them to be consistent within their la errors. We 
then varied the outer radius, i.e., 3" - 8 ", 3" - 12", 3" - 15", 
and 3” — 20” (maximum outer radius avoiding the interaction 
region, see Figure 7), performed the same analysis, and found 
that all the temperatures are consistent within their la errors 
(although in some cases the temperature in the half annular re¬ 
gion near the interaction region was poorly constrained). Two 
possibilities can explain this: (i) the shock has not reached the 
central ~10” region of G036N; (ii) the shock has passed the 
central ~10” region of G036N, but the core is dense enough not 
to be destroyed. We now show that possibility (ii) is less likely. 
We first applied the same method to G036S, but the data qual¬ 
ity is poorer and all the temperatures are consistent given the 
large uncertainties. Notice thatM 2 oo, >~ 200 , , c SO imd, and t r for 

G036N and G036S are all similar (Tables 3 and 4), so we expect 
that the shock travels similar distances in G036N and G036S at 
any given time. As a result, if the shock has reached the core 
of G036N, it should also have reached the core of G036S. The 
pressure ratio, p n /ps = 0.61 — 1.16, where p n and p s are the 
mean pressures in the central 3” — 10” region of G036N and 
in the central 0— 15" region of G036S respectively, although 
with large uncertainty, is consistent with 1. This implies that, 
if possibility (ii) is correct, the core of G036S should also have 
survived the merger induced shock, in contradiction with obser¬ 
vations. Thus, the consistent temperatures in the two half an¬ 
nular regions and in their parent regions in the core of G036N 
imply that the shock probably has not reached the core, so the 
lack of gas cooling in the core of G036S is unlikely caused by 
the merger disrupting a pre-existing cool-core in G036S. This is 
also consistent with our previous expectation of an early stage 
of the merger, suggesting that the simplified model in Section 
5.1 overestimates the time since the merger began. The actual 
time since the merger began is probably less than ~ 0.4 — 0.7 
Gyr (see Table 4). 

For cool-core clusters, some heating mechanisms must be 
taking place to prevent the gas from catactrophic cooling (see 
Fabian 1994 and references therein). Interestingly, the BCG of 
G036N hosts a radio source (Section 3.1), with a flux density 
of 19.7 mJy and a 1.4 GHz luminosity {L\ .4 ghz) of 1.32 x 10 24 
W Hz -1 after the X-correction (Singal et al. 2011) assuming 


The spectral index a is defined as S v v a , where S v is the flux density at frequency v. 
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a radio spectral index of 0.10 (Condon et al. 1998). The ex¬ 
tent of the radio source is comparable to the angular resolution 
of the NVSS (FWHM = 45", -120 kpc at the cluster redshift, 
see right panel of Figure 3). If it is an extended source, i.e., 
with radio lobes filled with relativistic particles and magnetic 
fields (e.g., Birzan et al. 2008), there may be undetected cavi¬ 
ties, which are generally thought to play a key role in the “radio 
mode” feedback (e.g.. Bower et al. 2006; Croton et al. 2006; 
Sijacki et al. 2007) in clusters of galaxies (see McNamara & 
Nulsen 2007, 2012; Fabian 2012; and references therein). Us¬ 
ing the scaling relation in Cavagnolo et al. (2010), we estimated 
the cavity power, P ca v ~ (0.8 — 1.4) x 10 43 erg s -1 , based on the 
1.4 GHz luminosity. The X-ray luminosity of the cool-core is 
(3.3 —4.9) x 10 43 erg s _1 , 2 — 6 times larger than P cav . Stud¬ 
ies of large samples show that although in some cases the cav¬ 
ity power is smaller than the X-ray luminosity, on average it 
is energetic enough to balance the cooling (e.g., Rafferty et al. 
2006; Birzan et al. 2008; Dunn & Fabian 2008; Cavagnolo et 
al. 2010). Higher resolution radio data are required to con¬ 
firm the extended nature of the radio source, and deeper Chan¬ 
dra observations are also necessary to identify any cavities in 
the core of G036N. Sun (2009) studied a sample of 161 BCGs 
and 74 strong radio AGNs (Ty 4 ghz > 10 24 W Hz -1 ) in 152 
nearby (z < 0.11) clusters and groups using Chandra archive 
data. He found that all BCGs with L \ .4 ghz > 2 x 10 23 W Hz -1 
in the sample have cool-cores. Although with a higher redshift, 
G036N is also consistent with his finding. 

To conclude, G036N hosts a small, moderate cool-core, 
while G036S has at most a very weak cool-core. This differ¬ 
ence is unlikely to be caused by the ongoing merger. G036N 
also hosts a central radio source, which may be heating the gas 
if the radio source is extended. Examination of the temperature 
variations in the core of G036N also suggests that the simpli¬ 
fied merger dynamical model presented in Section 5.1 overes¬ 
timates the time since the merger began, which should be less 
than ~ 0.4 — 0.7 Gyr. 

6. SUMMARY 

We present XMM-Newton and Chandra observations of 
PLCK G036.7+14.9 from the Chandra-Planck Legacy Pro¬ 
gram. The X-ray images reveal two close (in projection), yet 
clearly separated subclusters, G036N and G036S (Figure 1); 
spectral analysis yields similar redshifts for the two subclusters 
(Section 4.3) and a higher temperature between them (Section 
4.7), confirming that they are interacting. Excluding the inter¬ 
action region, the morphologies of G036N and G036S do not 
deviate significantly from spherical symmetry (Figure 1) and 
the surface brightness profiles can be well described by single 
p models beyond the central ~10" regions (Figure 5 and Ta¬ 
ble 1), suggesting that the merger is still at an early stage. The 
BCGs of the two subclusters are very close to the X-ray centers 
(Figure 3), consistent with an early stage of the merger, and in 
favor of a merger mostly along the line-of-sight. A simplified 
dynamical model suggests that the merger began about 0.8 Gyr 
ago, with an angle between the merger axis and the plane of the 
sky of ~ 80° (Section 5.1), in accordance with our expectations. 

Surface brightness profiles (Figure 5) and temperature pro¬ 
files (Figure 6) indicate that the gas may be cooling in the 
cores of both subclusters. By calculating the cooling time, we 
conclude that G036N hosts a small, moderate cool-core, while 
G036S has at most a very weak cool-core. Based on the tem¬ 
perature variations in the core of G036N, we suggest that the 


difference in core cooling times between G036N and G036S is 
unlikely to be caused by the ongoing merger (Section 5.3). The 
temperature variations also indicate that the simplified model 
(Section 5.1) overestimates the time since the merger began, 
which is probably less than ~ 0.4 — 0.7 Gyr. Interestingly, the 
BCG of G036N hosts an unresolved radio source (Figure 3). 
If the radio source is extended, the “radio mode” AGN feed¬ 
back may be taking place in the core of G036N. Using the 1.4 
GHz luminosity, we estimate a cavity power about 2 — 6 times 
smaller than the X-ray luminosity of the core of G036N. Higher 
resolution radio data and deeper X-ray data are necessary to 
study the nature of the radio source and the physical processes 
in the core, e.g., jet/ICM interaction. 

Since the merger is probably at an early stage, we apply the 
hydrostatic equation to the ICM in the regions opposite to the 
interaction region to obtain the total mass (Section 4.5). We 
treat the uncertainties in Hydrogen column density (Section 
4.2), redshift (Section 4.3), Chandra background subtraction 
methods (Section 4.1.2), and other XMM-Newton and Chandra 
cross-calibration uncertainties as systematic uncertainties (Sec¬ 
tion 4.6). We find that Chandra gives a slightly higher total 
mass than XMM-Newton , although the difference is not signifi¬ 
cant (Section 4.6). We also compare the X-ray mass to the SZ 
mass (Section 5.2). The SZ mass is higher than the X-ray mass 
of either subcluster, but is lower than the X-ray mass of the 
whole cluster, which is caused by the fact that Planck does not 
resolve PLCK G036.7+14.9 into two subclusters and interprets 
it as a single cluster. 

Despite the short exposure times of the Chandra and XMM- 
Newton observations and large systematic uncertainties, we 
were able to infer the basic merger kinematics and some use¬ 
ful time scales based on the morphologies and core properties, 
from which we further estimate the total masses of the con¬ 
stituent subclusters and the whole cluster. They are all self- 
consistent. Obviously, deeper X-ray data and multiwavelength 
observations as well as targeted simulations are required to 
study this interesting system in greater detail. 

As the angular resolution of Planck is substantially coarser 
than those of Chandra and XMM-Newton, it cannot resolve 
closely interacting subclusters or close cluster pairs. Because 
of the size-flux degeneracy, Planck has to use the X-ray de¬ 
termined position and radius to refine the SZ flux measure¬ 
ment (e.g., Planck Collaboration XXIX 2014). Much of this 
information comes from MCXC, which collects cluster param¬ 
eters from publicly available ROSAT All-Sky Survey based and 
serendipitous based catalogues (Piffaretti et al. 2011). These 
catalogues usually consist of observations too shallow to re¬ 
solve clusters in much detail. Consequently, clusters with small 
projected separations might not be resolved and are treated as 
a single cluster, which is then adopted by MCXC and Planck. 
If such close clusters occupy a considerable fraction, the con¬ 
sequence is that the Planck all-sky survey will overestimate the 
number of massive clusters. Using this biased mass function to 
constrain cosmological parameter will result in lower fly and 
higher ds. Thus, higher resolution X-ray observations by Chan¬ 
dra and XMM-Newton of the Planck cluster sample are neces¬ 
sary to identify the fraction of previously unresolved clusters 
like PLCK G036.7+14.9 and correct for this bias before ap¬ 
plying the Planck SZ mass function to constrain cosmological 
parameters. 
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Table 1 

Best Fitting /3 Model Parameters of the Surface Brightness Profiles 


Region 

Facility 

11 

r c (arcsec) 

r c (kpc) 

X 2 /d.o.f 

G036N 

XMM-Newton 

0.6418 ±0.0060 

14.57 ±1.36 

39.08 ±3.65 

118.8/94 


Chandra 

0.6635 ±0.0095 

14.65 ±1.22 

39.29 ±3.27 

128.6/54 

G036S 

XMM-Newton 

0.6850 ±0.0078 

23.65 ±1.34 

63.43 ±3.59 

121.3/94 


Chandra 

0.7104±0.0134 

22.43 ±1.66 

60.16 ±4.45 

123.4/44 


Table 2 

Best Fitting Hydrogen Column Density and Redshift 


Method 

Parameter" 

G036N 

G036S 

XMM-Newton 

N H b 

z 

0.136 (0.092™) 
o 141 5+0.0025 /n 14Q0+O.OO39\ 

0 1 5 3 5 +0 0058 (0 1 5 54 +0 0043 ') 

u.ijoj _ 00051 v^.ujh-_ 00059 ; 

Chandra baseline + 

back*0.95 + back*1.05 c 

N H b 

z 

0.136 (0.1641® 

0 1590+ 0 0057 in i <co4+0.005(X 

u.ioyU -0 0034 tti. loy^t—Q oo^e 7 

0 1639 +00018 to 1640 +0 0035 l 

W.iuOV— 0 0 o29 (U.lOW - 0 0015 J 

Chandra with soft f/ 

N H b 

z 

0.136 (0.1951® 

0 1602 +00052 (0 1S69 +0 0042 ) 
u.iouz _ 0 0054 (u.ijo ^_ 0 0050 ; 

O 1 ^9oF0-0048 ( r \ i z:q o-|-0.0027\ 
U.lOZU 0 0047 (U.lOJQ 0 0049 ) 


^Obtained by global spectral fitting with the regions shown in Figure 4. 

^Hydrogen column density, in units of 10 22 cm -2 . The value without errors was obtained from the relation given in Willingale et al. (2013), while the values with 
errors are the best fitting values from our analysis. Nh for G036S is always linked to that for G036N (see Section 4.2). 

c Chandra baseline background model and varying Chandra baseline background model by ±5%. Varying the background level by ±5% only affects regions with 
faint source emission. Since the regions used for the global spectral fitting cover the bright part of the cluster, we only used the baseline background to obtain the best 
fitting Nh and z- In deriving the temperature profiles, we used these Nh and z for the three Chandra background models. 

d Chandra baseline background model with a soft-band adjustment. 
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Table 3 

Radius, Gas mass, Total mass, and Gas Mass Fraction 


Region 

Method 

Parameter" 

A = 2500 

A = 500 

A = 200 

G036N 

XMM-Newton 

r& 

M gA 

M a 

/a 

0.52±0.01(0.55 ±0.02) 

1.42 ±0.28(1.51 ±0.32) 

2.25 ±0.18(2.69 ±0.29) 
6.32± 1.41(5.63± 1.38) 

0.99 ±0.03(1.03 ±0.03) 
3.00 ±0.62(3.12 ±0.70) 
3.20 ±0.28(3.60 ±0.29) 

9.42 ±2.21(8.69 ±2.16) 

1.38 ±0.04(1.43 ±0.04) 
4.33 ±0.91 (4.49 ± 1.02) 
3.46±0.32(3.85 ±0.32) 
12.57 ±3.02(11.71 ±2.97) 


Chandra 

baseline* 

''A 

M g , A 

M a 

/a 

0.55 ±0.03(0.53 ±0.03) 
1.51 ±0.32(1.48 ±0.31) 
2.77±0.46(2.52±0.38) 
5.49± 1.29(5.92± 1.36) 

1.10 ±0.06(1.07 ±0.05) 
3.22 ±0.69(3.19 ±0.68) 
4.44 ± 0.67(4.08 ± 0.59) 
7.29± 1.78(7.83± 1.89) 

1.56 ± 0.08( 1.51 ±0.07) 
4.65 ± 0.99(4.60 ± 0.99) 
5.02±0.73(4.61 ±0.65) 
9.31 ±2.32(10.02±2.49) 


Chandra 
with soft" 

'A 

Mg, A 

m a 

/a 

0.54 ±0.04(0.54 ±0.03) 
1.45 ±0.33(1.52 ±0.35) 
2.59 ±0.65(2.63 ±0.41) 
5.71 ± 1.53(5.81 ±1.40) 

1.20 ±0.10(1.15 ±0.06) 
3.47 ±0.77(3.46±0.80) 
5.75 ±1.47(5.04 ±0.73) 

6.21 ±1.76(6.91 ±1.75) 

1.85 ±0.23(1.66 ±0.10) 
5.51 ± 1.33(5.09 ± 1.21) 

8.80 ±3.55(6.07 ±1.11) 

6.70 ±2.38(8.49 ±2.35) 


Chandra 

back*0.95 rf 

r\ 

M sA 

M a 

/a 

0.54 ±0.03(0.52 ±0.03) 
1.46 ±0.31(1.42 ±0.30) 

2.55 ±0.41(2.30 ±0.34) 
5.77 ± 1.37(6.22 ±1.43) 

1.06 ±0.06(1.03 ±0.05) 
3.09 ± 0.66(3.03 ± 0.65) 
3.96±0.64(3.57±0.53) 

7.89± 1.99(8.50±2.09) 

1.50 ±0.08(1.45 ±0.07) 
4.46 ±0.94(4.37 ±0.93) 
4.43 ±0.69(4.01 ±0.58) 

10.12 ±2.61(10.94 ±2.75) 


Chandra 
back* 1 .OS' 1 

r& 

M sA 

Ma 

/a 

0.52 ± 0.02(0.54 ± 0.03) 
1.43 ±0.31(1.49 ±0.35) 
2.36 ±0.27(2.56 ±0.44) 
6.06± 1.41(5.87± 1.47) 

1.10 ±0.05 (1.08 ±0.07) 
3.21 ±0.74(3.22±0.77) 
4.37 ±0.55(4.20±0.75) 
7.40± 1.81(7.69±2.08) 

1.63 ±0.09(1.55 ±0.09) 
4.88± 1.16(4.70± 1.12) 
5.74 ±0.98(4.92 ±0.81) 
8.60 ±2.28(9.59 ±2.62) 

G036S/ 

XMM-Newton 

r A 

MgA 

Ma 

/a 

0.52 ± 0.02(0.56 ± 0.02) 
1.68 ±0.33(1.79 ±0.31) 

2.36 ±0.23(2.95 ±0.32) 

7.10± 1.42(6.07 ±1.13) 

0.98 ±0.03 (1.06 ±0.04) 

3.36 ±0.67(3.56 ±0.63) 

3.13 ±0.31(3.90 ±0.41) 
10.72±2.27(9.12± 1.77) 

1.37 ±0.05(1.48 ±0.06) 
4.74 ±0.94(5.02 ±0.90) 

3.38 ± 0.34(4.24 ± 0.47) 
14.03 ±3.04(11.85 ±2.39) 


a r^, M g a, Ma, and /a are the radius, gas mass, total mass, and gas mass fraction, measured inside a radius where the overdensity within that radius is A times the critical 
density of the Universe at the cluster redshift, while the units are Mpc, 10 13 M©, 10 14 M©, and %, respectively. The values outside and inside of the parentheses were 
obtained with Nh fixed at the total and best-fitting values (Table 2), respectively. 
b Chandra baseline background model, i.e., the blank-sky background. 
c Chandra baseline background model with a soft-band adjustment. 
d Chandra baseline background model varied by +5%. 
e Chandra baseline background model varied by —5%. 

f As G036S is close to the chip boundary in the Chandra FOV, we can only measure the temperature profile up to ~0. 87 * 500 - In this relatively small region, Chandra 
gives an essentially isothermal temperature profile, unlike the obviously declining temperature profile revealed by XMM-Newton at large radii where Chandra do not 
cover. Thus we did not measure mass-related quantities for G036S with Chandra. 
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Table 4 

Characteristics of the Shock and the Merger 


Region 

Parameter" 

XMM-Newton 

Chandra baseline* 

Chandra back*0.95 c 

Chandra back*1.05 d 

G036N 

J( e 

^sound^ 

tr S 

1.16±0.21 (1.19±0.30) 
1246+^4 ( 1353 + 54 ) 
0.94±0.18 (0.87±0.22) 

2.01 ±0.46 (1.91 ±0.38) 
1323+7[ (1 267+7j) 
0.57±0.14 (0.61±0.13) 

2.01 ±0.45 (1.90±0.38) 
I 29 O +74 (1232+gg) 
0.57±0.14 (0.61 ±0.13) 

2.02±0.46 (1.91 ±0.38) 
1349+g® 3 (1295+gg) 

0.59 ±0.14 (0.61 ±0.13) 

G036S 

J4 e 

^sound^ 

tr S 

1.25±0.19 (1.28±0.28) 
1224+H (1318+44) 
0.88±0.14 (0.86±0.19) 

2.27±0.44 (2.15 ±0.36) 
1295+|! (1244+!|q) 

0.52 ±0.11 (0.57 ±0.10) 

2.26±0.43 (2.14±0.35) 

127 1+62 ( 12 1 7 ±53> 

0.51 ±0.10 (0.56±0.10) 

2.28±0.44 (2.17±0.36) 
1317+$ (1266+®) 

0.53 ±0.11 (0.57 ±0.10) 


t h 

'merge 

e‘ 

0.78 (0.79) 

76 (78) 

0.80 (0.79) 

80 (80) 

0.80 (0.80) 

79 (79) 

0.80 (0.80) 

81(81) 


"Values outside and inside of the parentheses were obtained with Nh fixed at the total and best fitting values (Table 2), respectively. 
b Chandra baseline background model. 
c Chandra baseline background model varied by +5%. 
d Chandra baseline background model varied by —5%. 

^Mach number obtained using the temperatures in Region 1 and Region 2 (Region 3) in Figure 7 as the temperatures of the postshock and preshock gas in G036N 
(G036S). 

f Sound speed, in units of km s -1 , with the temperature measured in the [0.15 — 0.75] 7*500 region (Table 5). 

^Time for the shock to propagate to the subcluster center in units of Gyr. For G036S, as Chandra does not cover large enough region to measure the mass, hence row, 
we assumed that the ratios of 7*200 for G036N and G036S are the same for XMM-Newton and Chandra and that the relative errors are the same for both subclusters for 
Chandra , to obtain 7*200 of G036S for Chandra. 

/j Time in Gyr since the merger began, obtained from the simplified model in Section 5.1. 

'Angle in degree between the merger axis (assuming zero impact parameter) and the plane of the sky obtained from the simplified model in Section 5.1. 


Table 5 

Temperature, y Xi500 , and T S z,50o/>x,500 


Region 

Parameter" 

XMM-Newton 

Chandra baseline* 

Chandra back*0.95* 

Chandra back*1.05 d 

G036N 

kT e 

7x,50(/ 

g 07 + 0.52 i "7 q^+u./in 
"• y, -0.44 

1.82 ±0.43 (2.13 ±0.54) 

6 73 + u,w (6 17+ u /4 ) 

°- '"—0.72 V 0 - 1 -0.69 1 

2.17±0.52 (2.02±0.48) 

6.40+!!;^ (5. 84 + 0 ( 54 ) 
2.06±0.50 (1.91 ±0.46) 

7 00+ l u/ 

/.UU -0 72 

2.26 ±0.58 

( 6 - 45 +^;) 

(2.11 ±0.54) 

G036S 

kT e 

lx,500 7 

5.76j+® ( 6 . 68 +g;®) 
1.97 ±0.41 (2.24 ±0.42) 

6 - 45 +q|! (5.95+0(5!) 

2.11 ±0.48 (1.95±0.43) 

ft 71 + 0.62 7Q+0.61\ 

°— -0.61 l 3 - lu -0.50> 

1.99 ±0.44 (1. 86 ± 0.41) 

6 67+° " 
D, °— 0.61 

2.15±0.50 

(6.I6+3J9) 

(2.02±0.45) 

Whole 

YSZ,500/Y X ,500 8 

1.11 ±0.35 (0.96±0.30) 

0.98±0.33 (1.06±0.35) 

1.04 ±0.34 (1.12 ± 0.37) 

0.96 ±0.33 

(1.02±0.35) 


"Values outside and inside of the parentheses were obtained with Nh fixed at the total and best fitting values (Table 2), respectively. 
b Chandra baseline background model. 
c Chandra baseline background model varied by +5%. 
d Chandra baseline background model varied by —5%. 

^Temperature in keV measured in the [0.15 — 0.75]t*5oo region. For G036N, we adopted 7*5 00 = 1.0 (1.1) Mpc for XMM-Newton {Chandra), while for G036S, 7*5 00 = 1.0 
(1.0) Mpc fox XMM-Newton {Chandra), for simplicity. 

^Product of temperature and gas mass, where temperature was measured in the [0.15 — 0.75] 7*500 region and gas mass is that within 7 * 500 - lx,500 is in 10 14 Mq keV. 
5 See Section 5.2 for the definition of this ratio. Tsz ,500 h as been scaled to a redshift of 0.1547. 
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50.0 40.0 18:04:30.0 20.0 10.0 50.0 40.0 1804:30.0 20.0 10.0 

Fig. 1.— X-ray images (2.2 Mpc x 2.1 Mpc) of PLCK G036.7+14.9. North is up and East is to the left. Left : XMM-Newton EPIC (quiescent particle and 
residual soft proton) background subtracted and exposure corrected image in the 0.7 — 2.0 keV band, smoothed with a 15 77 Gaussian kernel. Each pixel has a size of 
2( 7 5 x 2'.'5. The large scale X-ray emission is elongated in a northeast-southwest direction, suggesting that PLCK G036.7+14.9 may be undergoing a merger. Two 
subclusters, G036N in the north and G036S in the south, were also resolved, although the separation between them is quite small (~72 77 = 193 kpc in projection). 
Outside the interaction region between the two subclusters, the morphologies of G036N and G036S do not deviate significantly from spherical symmetry. Right : 
Chandra 0.7 — 7.0 keV background subtracted and exposure corrected image, smoothed with a 3( 7 94 Gaussian kernel. The image is not binned; i.e., 1 pixel has a 
size of 0( 7 492 x 0( 7 492. The large scale X-ray emission exhibits the same behavior as XMM-Newton revealed, with more details of the two subclusters resolved. A 
point source is visible at the center of G036N, which is an AGN. A bow shaped gap (with an angle of ~145°), is seen between G036N and G036S, another indication 
of interaction between G036N and G036S. 



Fig. 2.— Left : Color-color {g — r vs. r — i) diagram of non-stellar objects in the FOV (7(7 x 7(7) of MOSCA. The shaded circle indicates the region from which 
galaxies were selected to derive a linear fit to the color-magnitude relation (right panel) of the cluster galaxies. Right : Color-magnitude (g — i vs. i) diagram of 
non-stellar objects in the same FOV as in the left panel. The dashed line is the best fitting color-magnitude relation for galaxies within the circle shown in the left 
panel, while the scatter (0.072 mag, la) about the best fitting relation is indicated by the dotted lines. Galaxies within twice the measured scatter around the best 
fitting relation, down to Iab ~ 21 mag, are represented by red dots and are considered to be the red sequence of the cluster. The BCGs (marked in Figure 3, left 
panel) of the two subclusters are located on the red sequence and colored green in particular in both panels. 
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Fig. 3.— Left : Optical mosaic image (1.2 Mpc x 1.2 Mpc) with the XMM-Newton 0.7 — 2.0 keV contours overlaid. The two blue squares (5 " x 5") mark the 
locations of the two BCGs identified in Figure 2. As can be seen, the two BCGs are very close to the X-ray centroids (defined to be the X-ray centers) of G036N 
and G036S. Right: NVSS 1.4 GHz radio image, same size as the left panel. The magenta crosses mark the X-ray centers of G036N and G036S. While there is a 
radio source (name labeled) at the location of G036N, no radio emission was detected at the position of G036S. The radio source should be hosted by the BCG of 
G036N. The green circle shows the angular resolution of the NVSS (FWHM = 45", 121 kpc at the cluster redshift). The extent of the radio source is comparable to 
the angular resolution of the NVSS. 



Fig. 4.— Regions used to perform spatial and spectral analyses. For either G036N or G036S, the central small green circle (35" in radius) represents the core, 
while the outer green ellipse represents the bright part of the subcluster. A third region, the ellipse exclude the core, is also defined. These three regions are the 
same for both XMM-Newton and Chandra. Spectra from these regions were extracted and fit together to determine the global properties (see Sections 4.1 —4.3) of 
G036N and G036S, with both XMM-Newton and Chandra. The magenta half circles mark 7*500 for both G036N and G036S, and temperature and surface brightness 
profiles were obtained within these regions (note that the outmost radii for these profiles are different from 7 * 500 ; see Figures 5 and 6). Left: XMM-Newton regions. 
For both G036N and G036S, 7*500 =1.0 Mpc (see Table 3). The outmost blue dashed annulus (480 // — 720" in radii) is largely cluster emission free, and, together 
with the six regions defined above for both G036N and G036S, is used to determine the various components of the background (see Section 4.1.2 for details). Right: 
Chandra regions. For G036N, 7*500 = 11 Mpc, while 7*500 = 10 Mpc for G036S (see Tables 3 and 5). The blue dashed region near the CCD chip is the region used 
to determine the parameters of the soft-band adjustment (see Section 4.1.1 for details). 
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R (kpc) R (kpc) 



Fig. 5.— Surface brightness profiles extracted from the half circular regions in Figure 4. Red data points with errors are the measured values. The green dashed 
lines are the best fitting single p models outside the central ~10 // regions, where there is excess emission relative to the p models. The p models generally fit the 
data very well in these regions. The black dotted lines in the central ~10 // regions are the interpolations of the p models. The outer boundaries are shown by the 
largest values of the x-axis. Upper Left : Surface brightness profile of G036N obtained by XMM-Newton. Upper Right : Surface brightness profile of G036S obtained 
by XMM-Newton. Lower Left : Surface brightness profile of G036N obtained by Chandra. Lower Right: Surface brightness profile of G036S obtained by Chandra. 
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Fig. 6 .— Temperature profiles extracted from the half circular regions in Figure 4. XMM-Newton data are marked by red cross, Chandra baseline background 
model by green diamond, Chandra baseline background model varied by —5% by blue square, Chandra baseline background model varied by ±5% by magenta 
circle, and Chandra baseline background model with a soft-band adjustment by cyan triangle. The horizontal error bars show the radial bin size. To show the data 
more clearly, we offset the Chandra data points by ±26.8 and ±13.4 kpc for the four Chandra background models. Upper Left : Temperature profiles of G036N 
obtained with Nh fixed at the total value, 0.136 x 10 22 cm~ 2 . Upper Right : Temperature profiles of G036N obtained with Nh fixed at the best fitting values, 
0.092 x 10 22 cm -2 for XMM-Newton, 0.195 x 10 22 cm -2 for the Chandra baseline background model with a soft-band adjustment, and 0.164 x 10 22 cm~ 2 for the 
other three Chandra background models. Lower Left'. Temperature profiles of G036S obtained with Nh fixed at the total value. Lower Right : Temperature profiles 
of G036S obtained with Nh fixed at the best fitting values. 



Fig. 7.— Regions used to determine the Mach number of the merger shock. Region 1 is the postshock region, and the symmetric region with respect to the 
subcluster center (magenta cross), Region 2 (3), is assumed to be the preshock region for G036N (G036S). Left : Regions overlaid on the XMM-Newton image. The 
lower right scale bar is 200 kpc. Right : Regions overlaid on the Chandra image. This image has the same size as the left panel. 
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Fig. 8.— Regions used to determine the mass for both X-ray and SZ measurements. The green solid circle, with a radius of rsz ,500 = 1.12 Mpc, is the region to 
measure the SZ mass by Planck, while the two magenta dashed circles, each with a radius of rx 5 oo> are the regions to “measure” (note that we only measure the mass 
in half circle and double it to obtain the total mass) the X-ray mass by XMM-Newton or Chandra. The green solid circle is centered at the Planck position, while the 
magenta dashed circles are centered at the X-ray centroids for both G036N and G036S. Left: XMM-Newton regions. For both G036N and G036S, Of ,500 =1.0 Mpc 
(see Table 3). The SZ region is larger by ~14% than the X-ray region (the union of G036N and G036S). Right : Chandra regions. For G036N, rx ,500 = 1.1 Mpc, 
while for G036S, rx 500 = 1.0 Mpc (see Tables 3 and 5 for details). The SZ region is larger by ~3% than the X-ray region. 
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APPENDIX 

TEMPERATURE, GAS MASS, TOTAL MASS, AND GAS MASS FRACTION PROFILES 

We present the model 2D and 3D temperature profiles, gas mass and total mass profiles, and gas mass fraction profiles in this 
Appendix. Sections 4.4 and 4.5 provide the details on how to obtain these profiles. 


Temperature Profile 



R(kpc) 



Gas Mass Faction Profile 



Fig. A9.— Temperature, gas mass and total mass, and gas mass fraction profiles for G036N obtained by XMM-Newton, with Nr — 0.136 x 10 22 cm -2 (total value). 
The vertical green dashed line marks 7*500 f° r ah three panels. Left : Temperature profiles. The blue solid and blue dotted lines are the best fitting 2D temperature 
profile and the corresponding lcr errors, while the red dash-dotted and red short-dashed lines are the best fitting 3D temperature profiles and the corresponding 
errors. The data points with errors are the measured 2D temperatures in each bin. Middle : Gas mass and total mass profiles. The upper red dash-dotted and red 
short-dashed lines are the total mass profile and the corresponding errors, while the lower blue solid and blue dotted lines are the gas mass and the corresponding 
errors. Right : Gas mass fraction profile. The blue solid and blue dotted lines are the gas mass fraction profile and the corresponding errors. 
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Fig. A11.— Same as Figure A9, except for G036S. 
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Temperature Profile 


Mass Profile 


Gas Mass Faction Profile 





Fig. A15.— Same as Figure A9, except for Chandra baseline background model varying by —5%. 
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Temperature Profile 
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Mass Profile 
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Fig. A16.— Same as Figure A9, except for Chandra baseline background model with a soft band adjustment. 


Gas Mass Faction Profile 
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Fig. A17.— Same as Figure A9, except for Chandra baseline background model and Nh — 0.164 x 10 22 cm 2 (best fitting value). 


Temperature Profile 


Mass Profile 


Gas Mass Faction Profile 





Fig. A18.— Same as Figure A9, except for Chandra baseline background model varying by +5% and Nh — 0.164 x 10 22 cm 2 (best fitting value). 





Fig. A19.— Same as Figure A9, except for Chandra baseline background model varying by —5% and Nh — 0.164 x 10 22 cm 2 (best fitting value). 
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Fig. A20.— Same as Figure A9, except for Chandra baseline background model with a soft band adjustment and Nh — 0.195 x 10 22 cm 2 (best fitting value). 























